Valley-isospin dependence of the quantum Hall effect in a graphene p-n junction 
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We calculate the conductance G of a bipolar junction in a graphene nanoribbon, in the high- 
magnetic field regime where the Hall conductance in the p-doped and n-doped regions is 2e^/h. In 
the absence of intervalley scattering, the result G — (e^//i)(l — cos"l>) depends only on the angle 
$ between the valley isospins {— Bloch vectors representing the spinor of the valley polarization) 
at the two opposite edges. This plateau in the conductance versus Fermi energy is insensitive to 
electrostatic disorder, while it is destabilized by the dispersionless edge state which may exist at a 
zigzag boundary. A strain-induced vector potential shifts the conductance plateau up or down by 
rotating the valley isospin. 

PACS numbers: 73.43.-f, 73.21.Hb, 73.23.-b, 73.50.Jt 



Recent experimentsi'^'^ have succeeded in fabricating 
junctions between p-doped and n-doped graphene, and 
have begun to investigate the remarkable properties pre- 
dicted theoretically.i^'^'i The conductance G of a p-n 
junction measures the coupling of electron-like states 
from the conduction band to hole-like states from the 
valence band, which in graphene is unusually strong be- 
cause of the phenomenon of Klein tunneling.'*'^ 

In the zero-magnetic field regime of Huard et alfi 
this coupling depends on the length scales character- 
istic of the p-n interface. In the high-magnetic field 
regime of Williams, DiCarlo, and MarcuSf2i the p-n junc- 
tion has a quantized conductance, which has been ex- 
plained by Abanin and LevitOA*^ as the series conduc- 
tance Gscrios = GpGn/ (Gp -\- On) of the quantum Hall 
conductances Gp,Gn in the p-doped and n-doped re- 
gions (each an odd multiple of the conductance quan- 
tum Go = 2e^//i). (The p-n-p junction experiments of 
Ozyilmaz et al.'^ are also explained in terms of a series 
conductance.) 

These results apply if the system is sufficiently large 
that mesoscopic fluctuations in the conductance can be 
ignored, either as a consequence of self-averaging by time 
dependent electric fields or as a consequence of suppres- 
sion of phase coherence by inelastic scattering?^ In a suffi- 
ciently small system mesoscopic conductance fluctuations 
as a function of Fermi energy are expected to appear. In 
particular, in the quantum Hall effect regime, the conduc- 
tance of a p-n junction is expected to fluctuate around 
the series conductance Gsorics m a small conductor (nano- 
ribbon) at low temperatures. 

In this paper we show that a plateau in the conduc- 
tance versus Fermi energy survives in the case of fully 
phase coherent conduction without intervalley scatter- 
ing. When both p-doped and n-doped regions are on the 
lowest Hall plateau (Gp = G„ = Go), we find a plateau 
at 
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G= iGo(l-cos$) 



(1) 



with $ the angle between the valley isospins at the two 
edges of the nanoribbon. A random electrostatic poten- 



N 



FIG. 1: Schematic top view of a graphene nanoribbon con- 
taining an interface between an n-doped and p-doped region 
(left panel) and between a normal (N) and superconduct- 
ing (S) region (right panel). Electron-like and hole-like edge 
states in the lowest Landau level are indicated by solid and 
dashed lines, respectively, with arrows pointing in the direc- 
tion of propagation. The electron-like and hole-like valley- 
polarized edge states hybridize along the p-n or NS inter- 
face to form a valley-degenerate electron-hole state. The two- 
terminal conductance G = GoT^h is determined by the prob- 
ability Teh that an electron-like state is converted into a hole- 
like state at the opposite edge (with Go = 2e^//i in the p-n 
junction and Go = 4e^//i in the NS junction). In the absence 
of intervalley scattering, Teh = ^(1 ~ cos "!>), with "I> the angle 
between the valley isospins of the electron-like state at the 
two edges.— 



tial is not effective at producing mesoscopic conductance 
fluctuations, provided that it varies slowly on the scale of 
the lattice constant — so that it does not induce inter- 
valley scattering. The dispersionless edge state that may 
exist at a zigzag edge (and connects the two valleys at 
opposite edges) is an intrinsic source of intervalley scat- 
tering when the edge crosses the p-n interface. The angle 
$ that determines the conductance plateau can be varied 
by straining the carbon lattice, either systematically to 
shift the plateau up or down, or randomly to produce a 
bimodal statistical distribution of the conductance in an 
armchair nanoribbon. 

Our analysis was inspired by an analogy between edge 
channel transport of Dirac fermions along a p-n interface^ 
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FIG. 2: Potential step at the p-n interface (with the shaded 
area indicating the energy range in the valence band). The 
electrostatic potential U{x) increases from to Uoo over a 
distance L around a; = 0. The Fermi level at Ef £ (0, Uoo) 
lies in the conduction band for negative x (n-doped region) 
and in the valence band for positive x (p-doped region). 



and along a normal-superconducting (NS) interface.— The 
analogy, explained in Fig. [1] is instructive, but it is only 
a partial analogy — as we will see. We present analytical 
results, obtained from the Dirac equation, as well as nu- 
merical results, obtained from a tight-binding model on 
a honeycomb lattice. We start with the former. 

The Dirac equation for massless two-dimensional 
fermions reads 



To (g) [v{p + eA) ■ a- + C/]* = 



(2) 



where E is the energy, v the Fermi velocity, p — 
{h/i){d/dx,d/dy) the canonical momentum operator in 
the x-y plane of the graphene layer, U{x) the electrostatic 
potential step at the p-n interface (shown in Fig.[2|), and 
A the vector potential corresponding to a perpendicular 
magnetic field B. The Pauli matrices cFi and act on 
the sublattice and valley degree of freedom, respectively 
(with (Tq and tq representing the 2x2 unit matrix). 

The Dirac equation ^ is written in the valley-isotropic 
representation, in which the boundary condition for the 
wave function 'J at the edges of the nanoribbon (taken 
at y = 0, W) has the form^ 



* = (j^ • x) (g) (sin Oa^ + cos 9 a^)^, 



(3) 



parameterized by an angle 6 and by the three- 
dimensional unit vector u on the Bloch sphere. The vec- 
tor u is called the valley isospin because it represents the 
two-component spinor of the valley degree of freedom,- 
An armchair edge has ■ z = 0, 6 = tt/2 (modulo 
tt), while a zigzag edge has ■ z\ = 1, 6 ~ (modulo 
tt). Confinement by an infinite mass has ■ z\ — 1, 
Q = -K jl (modulo tt). Intermediate values oiv ■ z and 
are produced, for example, by a staggered edge potential 
(having a different value on the two sublattices).'°!ii If 
the edge is inhomogeneous, it is the value of v and Q 
in the vicinity of the p-n interface (within a magnetic 
length Ini = ^Jh/eB from x = Q) that matters for the 
conductance. 
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FIG. 3: Dispersion relation E^(q) according to Eq. (|4]) of 
edge states near the Dirac point in the n region (solid curves) 
and in the p region (dashed curves). The color of the curves 
indicates the valley polarization (blue: -\-v, red; —v)- The 
three panels correspond to three different boundary condi- 
tions, and illustrate the transition from an armchair edge 
(leftmost panel) to a zigzag edge (rightmost panel). 



The boundary condition ^ breaks the valley degen- 
eracy of quantum Hall edge states r^ii^ii^ with different 
dispersion relations E'^ [q) for the two eigenstates | ± u) 
oi u ■ T. (We use the Landau gauge in which A is par- 
allel to the boundary and vanishes at the boundary. In 
this gauge the canonical momentum Tiq parallel to the 
boundary is a good quantum number.) In the n region 
(where J7 = 0) the dispersion relation is determined by 
the equations^ 



fE+ (q) = tan(0/2) , fE- {q) - -cotan {6/2), 



fE{q) ^ 



H^-2/2{qlm) 



£2/2-1 



{qlm) 



, S = Elm/hv, 



(4a) 
(4b) 



with Ha{x) the Hermite function. The dispersion rela- 
tion in the p region is obtained by E^ (q) ^ E^ (q) + Uoo ■ 

The dispersion relation near the Dirac point {E = 0) 
is plotted in Fig. [3] for three values of 0. (It does not 
depend on u.) For any 9 ^ Q (modulo tt) there is a 
nonzero interval lS.Ep of Fermi energies in which just two 
edge channels of opposite valley isospin cross the Fermi 
level (dotted line), one electron- like edge channel from 
the n region (blue solid curve) and one hole-like edge 
channel from the p region (red dashed curve). The case 

= is special because of the dispersionless edge state 
which extends along a zigzag boundaryfi^ As 6* ^ the 
interval l^Ep shrinks to zero, and at = (modulo tt) 
the electron-like and hole-like edge channels in the lowest 
Landau level have identical valley isospins. It is here 
that the analogy with the problem of the NS junction^ 
stops, because in that problem the electron and hole edge 
channels at the Fermi level have opposite valley isospins 
irrespective of 9. 

The two valley-polarized edge channels from the n and 
p regions are coupled by the potential step at the p-n 




FIG. 4: Dispersion relation at the p-n interface, calculated nu- 
merically from the Dirac equation for a step function potential 
profile. Each Landau level has a twofold valley degeneracy. 



interface. Edge states along a potential step which is 
smooth on the scale of the lattice constant a are valley 
degenerateji^iii because an electrostatic potential in the 
Dirac equation does not couple the valleys. The disper- 
sion relation, for the case of an abrupt potential step 
(a ^ L <C Im), is plotted in Fig. |H (It is qualitatively 
similar for L ^ Im ) The Fermi level now intersects 
with a two-fold valley degenerate edge channel of mixed 
electron- hole character. 

The two-terminal conductance of the p-n junction is 
given byi G = GoTeh, in terms of the probability Teh 
that an electron incident in an electron-like edge chan- 
nel along the left edge is transmitted to a hole-like edge 
channel along the right edge. We now show that this 
probability takes on a universal form, dependent only on 
the valley isospins at the edge, in the absence of inter- 
valley scattering. The argument is analogous to that in 
the NS junction,^ and requires that the electron-like and 
hole-like edge channels at the same edge have opposite 
valley isospins (±i/l for the left edge and zLvn for the 
right edge)jii^ 

Since the unidirectional motion of the edge states pre- 
vents reflections, the total transmission matrix itotai — 
tRtpntL from one edge to the other edge is the product 
of three 2x2 unitary matrices: the transmission matrix 

from the left edge to the p-n interface, the transmis- 
sion matrix tp„ along the interface, and the transmission 
matrix tn from the p-n interface to the right edge. In 



the absence of intervalley scattering ip„ 
portional to the unit matrix, while 



To IS pro- 



tx 



\ +iyx) {+1^x1 +e'^'-\-iyx) {-1^x1 (5) 



(with X = R) is diagonal in the basis | ± vx) of eigen- 
states of fx - t. The phase shifts (f)pn, (t>x, 4>'x need not be 
determined. Evaluation of the transmission probability 

Teh = \{+l^L\ttot.x\-l^R)? (6) 

leads to the conductance ([T]) with cos ^ — h>L ■ ^r- 
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FIG. 5: Conductance of an armchair nanoribbon containing 
the potential step U{x) — |[tanh(2a;/I/) + l]Uoo, calculated 
numerically from the tight-binding model in a perpendicular 
magnetic field {Im ~ 5 a). The step height Uoa is varied from 
below Ef (unipolar regime) to above Ef (bipolar regime), 
at fixed Ef — hv/lm and L — 50 a. The solid curves are 
without disorder, while the dashed curves are for a random 
electrostatic potential landscape {Kq = 1, ^ = 10a). The 
number J\f of hexagons across the ribbon is 97 (red curves), 
98 (blue), and 99 (green). The dashed horizontal line marks 
the plateau at G = | x 2e^ /h. 
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FIG. 6: Same as Fig. (5] for the case of a zigzag nanoribbon 
(A/" — 114 for the green curves and 115 for the red curves). 



To test the robustness of the conductance plateau to 
a random electrostatic potential, we have performed nu- 
merical simulations. A random potential landscape is 
introduced in the same way as in Ref. [l9l . by randomly 
placing impurities at A'^i,„p sites Ri on a honeycomb lat- 
tice. Each impurity has a Gaussian potential profile 
C/iexp(— |r — i2ip/2^^) of range ^ and random height 
Ui G (—(5,(5). We take ^ equal to the mean separation d 
of the impurities and large compared to the lattice con- 
stant a. The strength of the resulting potential fluctua- 
tions SU (r) is quantified by the dimensionless correlator 
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(7) 
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where the sum runs over all iVtot lattice sites in a 
nanoribbon of area A. 

Results are shown in Figs. [5] and [6] for an armchair and 
zigzag nanoribbon, respectively. The angle <i> between 
the valley isospins at two opposite armchair edges de- 
pends on the number M of hexagons across the ribbon: 
$ TT if TV is a multiple of 3, |$| = 7r/3 if it is not.^" 
Fig. [5] indeed shows that the conductance as a function of 
Uoo — Ep switches from a plateau at the ^-independent 
Hall conductance Gq in the unipolar regime {Uoo < Ep) 
to a ^-dependent value given by Eq. ^ in the bipolar 
regime (Uoo > Ep). The plateau persists in the presence 
of a smooth random potential (compare solid and dashed 
curves in Fig. [5]). By reducing the potential range we 
found that the plateaus did not disappear until ^ ^ 3 a 
(not shown). 

As expected in view of the intervalley scattering pro- 
duced by the dispersionless edge state in a zigzag nano- 
ribbon, no such robust conductance plateau exists in this 
case (Fig.[S]). In the presence of disorder the conductance 
oscillates around its ensemble average Gq/^, in a sample 
specific manner. The numerics for any given realization 
of the disorder potential satisfies approximately the sum 
rule G{J\f) + ClAf + 1) Go, ioT which we have not yet 
found an analytical derivation. 

The valley-isospin dependence of the quantum Hall ef- 
fect in a p-n junction makes it possible to use strain 
as a means of variation of the height of the conduc- 
tance plateaus. Strain introduces a vector potential term 
evTz (8) {5A ■ (t)^' in the Dirac equation correspond- 
ing to a fictitious magnetic field of opposite sign in the 
two valleys f^i^iS^!^ This term rotates the Bloch vector 
of the valley isospin around the z-axis, which in the case 
of an armchair nanoribbon corresponds to a rotation of 
the valley isospin in the x-y plane. Strain may appear 
locally at an armchair edge by passivation of the carbon 



bonds. W (The resulting change St of the hopping energy 
T changes $ by an amounl^ 5$ = 2^/3 St /t.) Random 
strain along the p-n interface, resulting from mesoscopic 
corrugation of the carbon monolayer,^** corresponds to 
a random value of the angle <f> in the conductance for- 
mula ^ . A uniform distribution of $ implies a bimodal 
statistical distribution of the conductance, 

P{G) = - f d$5[G- iGo(l-cos$)] 

= [tt^G{Gq - G)]-^/'\. 0<G<Go, (8) 

distinct from the uniform distribution expected for ran- 
dom edge channel mixingf^ 

In summary, we have presented analytical and numer- 
ical evidence for the existence of a valley-isospin depen- 
dent conductance plateau in a p-n junction in the quan- 
tum Hall effect regime. In recent experiments^'^ the 
conductance was simply the series conductance of the 
p-doped and n-doped regions, presumably because of lo- 
cal equilibration. We have shown that the mesoscopic 
fiuctuations, expected to appear in the phase coherent 
regime;^ are suppressed in the absence of intervalley scat- 
tering. The conductance plateau is then not given by the 
series conductance, but by Eq. ([1]). The same formula 
applies to the conductance of a normal-superconducting 
junction in graphene;^ revealing an intriguing analogy 
between Klein tunneling in p-n junctions and Andreev 
refiection at NS interfaces.— 121 
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